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Standard approximations for the exchange-correlation functional have been found to give big errors 
for the linearity condition of fractional charges, leading to delocalization error, and the constancy 
condition of fractional spins, leading to static correlation error. These two conditions are now unified 
for states with both fractional charge and fractional spin: the exact energy functional is a plane, 
linear along the fractional charge coordinate and constant along the fractional spin coordinate with 
a line of discontinuity at the integer. This sheds light on the nature of the derivative discontinuity 
and calls for explicitly discontinuous functionals of the density or orbitals that go beyond currently 
used smooth approximations. This is key for the application of DFT to strongly correlated systems. 
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The great success of density functional theory (DFT) 
[T] is clouded by spectacular failures ^ that can mani- 
fest themselves in a variety of problems, from molecular 
bond stretching to magnetism and the band-gap of ma- 
terials. Some of these failures in practical calculations 
are even considered a breakdown of the theory itself, as 
exemplified by the concept of strong-correlation, which 
often refers to the supposed qualitative collapse of the 
single particle picture and hence DFT. In this Letter we 
present an important condition of the exact DFT which 
is missing in currently used approximations and is vital 
to understand strongly-correlated systems, in particular 
their bandgaps. 

The fundamental gap of a N-electron system with ex- 
ternal potential v{r) is given by the difference of the ion- 
ization energy, /, and electron affinity. A, 

4ap^°'' = {Ev{N - 1) - E.,{N)}^{E,,{N) ~ E.,{N + 1)} 

(1) 

or, due to the straight-line behavior of the exact func- 
tional [3; , by a difference of derivatives at iV 
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Consider a prototypical example of a strongly corre- 
lated system, a stretched H2 molecule. At its dissociation 
limit, there are two fractional spin hydrogen atoms each 
with half a spin up [a) electron and half a spin down {(3) 
electron, H[i, i]. This is a one-electron system that is 
degenerate in energy |5j with the normal hydrogen atom 
with one a electron, H[1,0], and therefore has the same 
£;i^tpCS°'-(= QA72Eh). However H[i, \] has identical a and 
(3 orbitals by symmetry, and hence all available function- 
als give i?g°p^ = 0. Where is the gap? What is missing 
from all available exchange-correlation functionals? In 
other words, what is the nature of the exact derivative 
discontinuity ^5, 6J? 

In this Letter we combine fractional charges and frac- 
tional spins to give a unified condition for the exact func- 
tional. This clarifies the behavior of the energy as the 



number of electrons passes through an integer and high- 
lights an important discontinuity in the energy expres- 
sion. An illustrative functional for the hydrogen atom is 
developed which clearly shows this discontinuous behav- 
ior and moreover how it manifests itself, with a discon- 
tinuous derivative. 

We follow the methodology of Yang, Zhang and Ayers 
[7J and examine systems at their dissociation limit. First, 
we start with the simplest case. Consider the hydrogen 
molecular ion HJ, which has one electron, one proton 
at site Ri and another proton infinitely far away at R2. 
The one electron can be at either of the two sites, and 
it can be spin up (m^ = |) or spin down {mg — —\)- 
Thus there are four degenerate electronic ground states, 
^ms,i — '^ms {Ri) I where {Ri) is the ground state 
of a hydrogen atom located at position Ri with the spin 
projection wLs and energy E(H). Consider the equally- 
weighted wavefunction 

^ = Y. *™..'/^- (3) 

The electron density of this wavefunction is 
2 2 ^ 

where Pm^^i is the density of ^'m.^./. The density p of 
Eq. Q is simply the sum of two identical densities 

= J (^pij + p_i_i^ ,1 = 1,2, separated from each 

other by an infinite distance. The energy of ^ and p 
is E{H). While the density p is ti-representable, is 
non-u-representable [U [9], and is the density of an iso- 
lated subsystem inside the supermolecule described by 
the wavefunction of Eq. ^ . 

Now we consider the behavior of the exact energy func- 
tional, Ey [p] , for this density. E^ [p] possesses the follow- 
ing properties: (A)^.;; [p] is exact for any (pure-state) v- 
representable density. Hence, for the total density in Eq. 
([4]), we have Ey[p\ = E{H). {^)E^ [p] is size extensive. 
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Therefore, E^[p\ = X^Li -^"[p'l- (C)-Eu [p] is translation- 
ally invariant. Therefore, = Ey[p^]. From (A),(B) 
and (C) it follows that 

E,[p'] = \E{H). (5) 

Thus, the exact energy for the non-w-representable den- 
sity = \{pi^i + \^{^)- This density has a 
fractional charge of | and fractional spins of j up-spin 
and \ down-spin. Its energy is ^E{H), the average of 
hydrogen atoms with 1 and electrons, independent of 
the fractional spins. Thus in this specific case, we obtain 
a condition for the exact functional for a density with 
fractional charge and spins, which extends the previous 
results of fractional charge [3] and fractional spin |1] . 

We now generalize our result of Eq. (|5| to include 
general fractional charge and fractional spins. We also 
extend the discussion to general degeneracies, instead 
of just a degeneracy because of spin symmetry. Con- 
sider an external potential v{r) that has two sets of de- 
generate grounds states: A^-electron degenerate ground 
states with energy Ey{N), wavefunctions {^N,is,'i = 
1, 2, g^r) and densities {pN,i,'i = 1, 2, ...,g7v), and 
{N + l)-electron degenerate ground states with energy 
Ey{N + 1), wavefunctions {^N+ij,j — 1, 2, gAr+i) 
and densities {pN+i,j,j — 1, 2, (Jtv+i)- For the den- 
sity P = \ Y!iZiCiPN,i + \ Y.%V d.jPN+i,j where {a} 
and {dj} are positive and finite integers, and satisfy 
the normalization conditions q — J2i=i'^i + 12j=i^ 
p = d-j, and q — p — J2i=i the exact energy 

functional satisfies the following equation 
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S^E.,iN) + ^^E,iN+l). 



(6) 

Eq. ([5]) is a special case of the general result, Eq. (|6|, 
and the proof is given in [10]. Eq. ^ is also valid in 
first-order reduced density-matrix functional theory. 

We analyze the simple case of a hydrogen atom with 
general spin up and spin down occupations, H[nQ,n^], 
to illustrate the scenario of fractional charges and frac- 
tional spins combined. This is of key importance for the 
consideration of strongly correlated systems and more 
especially the band-gap of Mott insulators in the non- 
magnetic phase. We focus on the simple strongly cor- 
related system H[^, ^J, which can be viewed as the in- 
finitely stretched limit of H2 3] or as the infinitely ex- 
panded limit of a lattice of hydrogen atoms [TT] but with 
zero spin-density everywhere. We now address the very 
interesting question: Is there a gap for H[i, ^]? 

The answer to this problem in terms of the total en- 
ergy can be understood from Eq. Fig. 1 shows the 
behavior of the exact energy functional for the hydro- 
gen atom with < tIq, < 1, < n/3 < 1, and < n < 2 
{n — Ua+np). The exact functional has flat plane behav- 
ior, linear along the fractional charge coordinate and con- 
stant along the fractional spin coordinate. This gives two 
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Figure 1: Exact energy (a.u.) for densities with fractional 
charges and fractional spins in H[nQ,n^]. Fractional charges 
occur along the green lines and fractional spins arise at the 
intersection ric, + Ufs = 1. The line for Ua = n/3 is highlighted 
in blue. 



flat planes, one for < n < 1 and other for 1 < n < 2, 
that intersect with a discontinuity at n = 1. The sim- 
ple fractional charge states correspond to the edge lines 
connecting [0,0] with [1,0] or [0,1] and then with [1,1], 
and the question of the gap in H[i, |] is highlighted by 
the blue line connecting [0,0] with [|, ^] and then with 
[1, 1]. If we now analyze this problem from a total en- 
ergy perspective it is clear that the energy of H[i, i] is 
degenerate with the normal H[1,0] atom, and also upon 
addition and removal of an electron (or any infinitesimal 
amount of an electron) the energy change is again exactly 
the same as the normal atom. This means that H[^, |] 
has an energy and derivatives, and therefore gap, 

that are exactly the same as H[1,0]. 

To gain more insight examine the behavior of sev- 
eral approximate functionals for the energy of H[nQ,,n^] 
shown in Fig. 2. If we first consider the fractional charge 
behavior, Hartree-Fock (HF) shows the exact straight 
line behavior between [0,0] and [1,0] as HF is exact for 
one orbital systems. However, it exhibits an incorrect gap 
for H[l, 0] due to concave curvature for 1 < n < 2, fur- 
ther characterized as localization error [12j . and incorrect 
energy of H[l, 1]. The behavior of semilocal functionals 
is exemplified by BLYP in Fig. 2b. They have an incor- 
rect convex behavior for fractional charges which usually 
leads to the underestimation of / and the overestimation 
of A from the derivatives [T3] and a delocalization error in 
larger systems [H]. The MCY3 functional has been de- 
veloped to give improved behavior for fractional charges 
[14] . so it greatly improves the edge lines in Fig. 2c and 
hence the gap from the derivatives for H[l, 0]. For all of 
the above functionals there is a failure to correctly predict 
the energy of H[i, ^] associated with the static correla- 
tion error of fractional spins However much more 
interesting is the gap of this system. This is most clearly 
illustrated in Fig. 2e where only the energy of the line 
"-a — is shown for < n < 2. H[i, i] is in the middle 
at n = 1 and the derivative gap is given by the difference 
between the slope to the left and the slope to the right. 
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Figure 2: The same as Fig. 1 for approximate functionals. All calculations are self-consistent using a cc-pVQZ basis set. 



It is clear that HF, BLYP and MCY3 all have a massive 
failure as they incorrectly have no gap for this system. 
This is a direct consequence of the smooth continuous 
behavior of all these different exchange-correlation func- 
tionals in terms of the orbitals. 

To understand the discontinuous behavior that is 



needed to give the correct gap, a simple illustrative func- 
tional is built to attempt to describe the whole two plane 
behavior of H[nQ, np] using a linear exchange term, a lin- 
ear correlation term and the challenging linear correction 
of the quadratic Coulomb term; 



1aEx[pa] + fcEc[Pa, Pjs] + (nhigh " n)J[p'^'' 



|-|-(n-niow)J[pa' 



"high] 

'Pf3 J 



'J[naPa,nf3Pf3], (7) 



with the general repulsion J[pa,Pb] = / /(/Oa(r) + 
P6(r)(pa(r') + p{,(r'))/|r — r'jdrdr'. There is an exphcit 
discontinuity in the correlation and Coulomb correction 
terms such that if ?^ < 1 /c = 0; njow = 0, p2'°'" = 
0,nhigh = l,/?"'""'' = riaPa/n and if n > 1/c = 
(n - l);niow = l,/?^'"" = (1 ^ n„,)p„/{2 - n),nhigii = 
2,^;^""=" ^ p„. We use HF for E^[p„] and LYP for 
Ec[pa,pf}\- Note that irrespective of the form used for 
Ec it gives exactly zero correlation energy for any one- 
electron system due to the prefactor, fc- Eq. ([T]) is a 
functional of the orbitals and occupation numbers and 
can be viewed within reduced-density-matrix functional 
theory, which has also been used to tackle strongly cor- 



related systems )15J . but in this work we only consider 
optimization of the orbitals as in standard DFT. This 
idea can be generalized and placed much more clearly 
within DFT if the normalization of the exact-exchange 
hole [ini E] is used instead of the occupation number. 
The performance of this functional in Fig. 2d shows 
a qualitative improvement over normal functionals and 
now resembles the exact functional of Fig. 1. It recovers 
the overall feature of two intersecting planes with a com- 
pletely flat plane between < n < 1 and an almost flat 
plane between 1 < n < 2, though slightly curved due to 
the approximate nature of the dynamic correlation. 

To examine the gap of Eq. (|2]) , consider the derivatives 
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Table L Energy gap (a.u.) from the derivatives (Eq. [sjl. 



BLYP HE MCY3 Eq. (7l Expt 



H[1,0] 



I- A 0.250 
I 

A 
I- A 
I 
A 



0.546 0.425 
0.272 0.500 0.448 
0.022 -0.046 0.023 

0.252 
0.252 



0.576 
0.500 



0.472 
0.500 





0.239 0.227 
0.239 0.227 



-0.076 0.028 

0.576 0.472 

0.500 0.500 

-0.076 0.028 



of the energy expression, 
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dN 
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dN 
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Ae 



A. 



-GKS 



-V^c (8) 



where Ae^^^ is the difference in the generalized Kohn- 
Sham frontier eigenvalues and incorporates all the dis- 
continuity due to a change of orbitals T3], including the 
smooth part, A^c, of an orbital dependent exchange- 
correlation term. Furthermore V^c represents only the 
explicit discontinuity of the exchange-correlation term, 
and hence goes beyond previous work [131 EH]- To un- 
derstand the nature of the exact derivative discontinuity 
[51 |S] both Aa-c and have to be considered. The HF, 
BLYP, MCY3 and all known approximations are only 
smooth functionals of the orbitals and therefore Vxc is 
zero. Only some methods based on the Hubbard model 
have a non-zero P^jc [19]. The new functional, Eq. Q 
has an explicit derivative discontinuity at n = 1 

V^^-^ = ~2J[pIpI] + JIpIp}] + E^^^IpIp}]. (9) 

Table 1 shows /, A and the gap from the derivatives. 
The gap of H[l, 0] is underpredicted by BLYP due to its 
convex behavior and overpredicted by HF due to its con- 
cave behavior and bad endpoint, while MCY3 improves 
in this case due to its straight line behavior. However 
none of these normal functionals offer a gap for H[i, ^] 
due to a lack ofV^c- The functional in Eq. ([t]) reveals the 
correct picture with a reasonable gap that is the same for 
both H[i, i] and H[1,0], although A is slightly underpre- 
dicted due to the concave nature for 1 < n < 2. Moreover 
these two gaps are both due to the explicit discontinu- 
ity of the functional, V^c, and there is no contribution 
from the orbital discontinuity, as the a and P orbitals 
from the self-consistent solution of Eq. ([t]) are always 
degenerate at n = 1. This raises an interesting ques- 
tion about normal functionals which erroneously break 
this degeneracy for n = I when Ua =/= They give a 
gap that is therefore of the wrong nature, due to a dis- 
continuity because of the orbital dependence, which only 
mimics the correct explicit discontinuity with degenerate 
orbitals. This erroneous behavior of normal functionals 
must no doubt contribute to the incorrect prediction of 



quantities related to these orbitals such as hyperfine, spin 
and magnetic properties. 

In this Letter a new exact condition for the energy 
functional is derived that shows the combined behavior 
for fractional charge and spins. It shows a discontinuous 
behavior when passing through the integer that reveals 
the explicitly discontinuous nature of the derivative of 
the exchange-correlation functional. This is most clearly 
highhghted by the gap of H[i, ^]. A gap for this system 
only appears if the exchange-correlation functional has 
an explicit discontinuity, because the orbitals are degen- 
erate and give no contribution. For Mott insulators, the 
unit cell has fractional spins corresponding to different 
magnetic phases and the correct gap prediction critically 
depends on this explicit discontinuity. For the future it 
is crucial to develop better approximations that go be- 
yond smooth orbital functionals by including an explicit 
discontinuity, along the line suggested by Eq. ([t]). This 
is essential for the advancement of DFT towards the cal- 
culation of strongly correlated systems. 
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APPENDIX: PROOF OF EQ. (6) 

We follow the methodology of Yang, Zhang and Ayers 
[7!| . Consider a supramolecular system with the following 
external potential Wtotai(r) = J21=i ^(i" ~ R-Oi namely it 
has q copies of the potential t;(r) each located at a site R; 
which are infinitely far away from each other. There is a 
total of qN + p electrons {N, p and q are all positive and 
finite integers and q > p). Since the sites are separated 
by infinite distances, the total system is simply composed 
of q subsystems in identical external potentials v{r), with 
no interaction between the subsystems. Its ground state 
has (q—p) A^-clectron subsystems and p (A^ + l)-electron 
subsystems, assuming the convexity condition 

E,{N) <{E,{N+l)+E,{N-l))/2, (10) 

which is known to hold for atoms and molecules from 
experimental data [3J|20J. The ground state is degenerate 
and its energy is 



iq-p)E,iN)+pE,iN + l) 



(11) 



The total ground-state wavefunction is an antisymmet- 
ric product of q separated ground state wavefunctions. 
One possible state is following: The first p locations, 
Ri...Rp, each have {N + 1) electrons; within these p 
locations, the first di sites have the degenerate wave- 
function <I>Ar+i_i, the second ^2 sites have the degenerate 
wavefunction ^n+i,2, and the last dg^^-^^ sites have 
the degenerate wavefunction ^N+i,gN+i- In this way, 
p = J2^jZV '^j- The remaining q—p locations, Rp+i...Rg, 
each have N electrons; within these q — p locations, the 
first ci sites have the degenerate wavefunction ^n.i, the 
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second C2 sites have the degenerate wavefunction ^n,2, ■■■ 
and the last Cg„ sites have the degenerate wavefunction 
^N.gN- III this way, q~p — c-i- Then this state has 

the wave function 



^-1 = A{$Ar+i,i (Ri) ...4>Ar+i,i (RdJ 

*i'Af+l,2 (R-di + l) (R-di+da) 

'^NS (R-p+l) ■■■'^NA (R-p+ci) 

'^Na (R-p+ci + l) (R-p+ci+cs) 



■} (12) 



Permutation of any two locations with different states 
{'^NA, or <i>Ar+ij) generates a different {qN + p)- 
electron wavefunction. There are a total of m = 
'''/(Hi" ciiJI'^'+i dji) such degenerate wave functions. 

For any wave function ^'fe, a particular site R/ can ei- 
ther have the wavefunction '^N,i, or <i>7v-|_ij. In all such 
wavefunctions {^k, k — 1, .., to}, the number of times any 
location Rs having the wave function <i>jv+i.n is equal to 

the corresponding number for ^N,n is equal to TOjv,n = 
('?-i)!/((d„-i)!f]J'" ci!f]^'^+i dj!) = ""irx/g. In analogy to 
Eq. (3) of the paper, the following equally-weighted 
wavefunction is also a degenerate wavefunction 

VTO — ' 

the density of which is 

9 I I 9N+1 \ 

1=1 ^ \i:=i j=i J 

In this particular state, with the degenerate energy of Eq. 



(Ill, all the q subsystems have the same electron density 



except by translation. Following the arguments leading 
to (A), (B) and (C) of the paper we have 



y--PE,{N) + ^E.,{N + l) 
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